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Long-Range Atom—Wall Interactions and Mixing Terms: Metastable Hydrogen 
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We investigate the interaction of metastable 2 S hydrogen atoms with a perfectly conducting wall, 
including parity-breaking S—P mixing terms (with full account of retardation). The neighboring 
2 Pi /2 and 2 -P 3/2 levels are found to have a profound effect on the transition from the short-range, 
nonrelativistic regime, to the retarded form of the Casimir-Polder interaction. The corresponding P 
state admixtures to the metastable 2 S state are calculated. We find the long-range asymptotics of 
the retarded Casimir-Polder potentials and mixing amplitudes, for general excited states, including 
a fully quantum electrodynamic treatment of the dipole-quadrupole mixing term. The decay width 
of the metastable 2 S state is roughly doubled even at a comparatively large distance of 918 atom 
units (Bohr radii) from the perfect conductor. The magnitude of the calculated effects is compared 
to the unexplained Sokolov effect. 
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Introduction. — The investigation of atom-wall interac¬ 
tions for atoms in contact with conducting materials has 
a long history. Starting from the works of Lennard- 
JonesQ, Bardeen Q, Casimir and Polder Q, and Lif- 
shitz Q, research on related matters has found contin¬ 
uously growing interest over the last decades @-@|- In 
the non-retarded regime (close range), the interaction en¬ 
ergy scales as 1/Z 3 with the atom-wall distance Z , while 
for atom-wall distances large in comparison to a typi¬ 
cal atomic wavelength, the interaction energy scales as 
1/Z 4 (see Chap. 8 of Ref. 0)- The leading term is given 
by virtual dipole transitions, while multipole corrections 
have recently been analyzed in Ref. 0- The symmetry 
breaking induced by the wall leads to dipole-quadrupole 
mixing terms, which lead to admixtures to metastable 
levels [Tl], 0| . While this effect has been analyzed in 
the non-retarded van-der-Waals regime 0 , 0 ] , a fully 
quantum electrodynamic calculation of this effect would 
be of obvious interest. 

This fact is emphasized by the curious observa¬ 
tion of a long-range, and conceivably super-long-range 
(micrometer-scale) interaction of metastable hydrogen 
2 S atoms with a conducting surface (the so-called 
Sokolov effect, see Refs. jl3Hl6| ). It is not far-fetched 
to suspect that this effect could be due to a quantum 
electrodynamically induced tail of the dipole-quadrupole 
mixing term in the atom-wall interaction. Namely, for 
the hydrogen 2 S atom, the neighboring 2P 1 / 2 and 2P 3 / 2 
levels are removed only by the Lamb shift and fine- 
structure, respectively, while it is known that virtual 
states of lower energy can induce long-range tails in atom- 
wall interactions, as well as in the Lamb shift between 
plates (see Refs. 0 - 0 ]). The large admixtures typi¬ 
cally induced in atomic systems when a metastable level 
couples to nearly degenerate states of opposite parity sug¬ 
gest that a closer investigation of the hydrogen system is 
warranted. Atomic units with h = 4-7reo = 1 and c = 1/a 
are used throughout this Rapid Communication, where 
a is the fine-structure constant. The electron charge is 
explicitly denoted as e unless stated otherwise. 

Retardation of the atom-wall interaction .—The quan¬ 


tum electrodynamic (QED) length-gauge interaction 

Hr = -e r ■ E - r 3 dE^dr 1 + ... , (1) 

follows naturally from the formalism of long-wavelength 
QED interaction Hamiltonian 0 , 0 ] (f denotes the elec¬ 
tron coordinate). In contrast to the vector potential, the 
electric field strength (operator) is gauge-invariant (this 
point has given rise to some discussion, see Ref. 0) and 
reads as [cf. Eq. (2.3) of Ref. 0 ] 

OO _ 
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where r = ry + ze z with ry = xe x + ye y , while fey = 
k x e x + k y e v , also k± = k z e z , and L = |fej_ |- The com¬ 
mutator relation is [a s (fcy,P), a)[,(fcy,P)] = S ss ' <5*- 2 )(fcy — 
fc') S(L — L') for the annihilation and creation operators 
a s and a\. In order to evaluate the interaction Hamilto¬ 
nian CD, one shifts z —> Z + z where Z is the coordinate 
of the atom’s center (nucleus). The proton is at (0,0, Z), 
while the atomic electron coordinates are ( x,y,Z + z). 
The surface of the perfect conductor is in the xy plane, 
i.e., in the plane described by the points (x,y, 0). The 
unperturbed Hamiltonian Hq is the sum of the free ra¬ 
diation held and the unperturbed atom [see Eq. (2.1) 
of Ref. 0 and Eq. (3.2) of Ref. 0]. For a refer¬ 
ence ground state |n), second-order perturbation theory 
leads to a known result given in Eq. (8.41) of Ref. [9] or 
Eq. (27) of Ref. 0, which involves the symmetric sum 
with imaginary frequency in the argument of the dynam¬ 
ics polarizability n(±iw). The Wick rotation of the vir¬ 
tual photon integration contour, leads to the symmetriza- 
tion i w <->■ —\oj but cannot be done for excited reference 
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states. We use second-order perturbation theory to eval- obtain [cf. the discussion following Eq. (2.12) of Ref. [20|]. 
uate A E = (n |(—er ■ E ) [1 /(E n — Hq)] (— er ■ E)\n) and 
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where the identity f R3 d 3 fc = 2 f 2w d<£ d L f^° dww , with w = + k 2 , and L = \k z \ has been used in order to 

transform the integration measure. The virtual states are denoted as Ig), and their energy difference to the reference 
state is denoted as £ q = E q — E n . In contrast to the velocity gauge [2C|, there is no seagull term to consider, and it 
is not necessary to add the electrostatic interaction with the mirror charges by hand }3l| . It is an in principle well 
known (see the remarks following Eq. (A.22) in Appendix A of Ref. (Hi), but sometimes forgotten wisdom that the 
Coulomb interaction does not need to be quantized in the velocity gauge [3l| • The integration with respect to ui leads 
to logarithmic terms [see the Appendix of Ref. [20j], After the subtraction of Z-independent terms (the subtraction 
is denoted by the = sign), one obtains 
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T(x) = sin(x) Ci(x) - cos(x) Si(x) + | cos(x), X = 2|£ q | Z , e(£ q ) = Q(£ q ) - ©(-£,) . (4c) 

Here, Ci(x) = - dt and Si(x) = f 0 X dt^M ; and U(x) = ^T( X ), while T( X ) = x _1 - -§^U{ X ). We confirm 
the result given in Eq. (2.18) of Ref. [2C| and represent the “distance-dependent Lamb shift” as 
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We should perhaps clarify that the Z-independent contribution to the Lamb shift (the ordinary “free-space Lamb 
shift”) is absorbed in the subtraction procedure denoted here by the “=” sign in Eqs. 0, ©, © and (18D . The 
Z-dependent position of the energy level is obtained after adding the “free-space Lamb shift” C and “free-space fine 
structure” T given in Eq. (fT^l) to the Z-dependent energy shifts given in Eqs. m and 0. In the nonretardation 
limit, the Z-dependent results given in Eqs. m and © are replaced by the respective terms of the nonretarded 
potential m- This (somewhat subtle) point is not fully discussed in previous works on the subject fl7l - [2l| and 
therefore should be mentioned for absolute clarity. 

The term —X -2 i n the coefficient multiplying |(n|rj||g)| — 2 |(n|z|(?)| 2 vanishes after summing over the entire 
spectrum of virtual states; it is obtained naturally in the length gauge and otherwise cancels a term in the expansion 
of the energy shift for large X (even before the application of the sum rule, which is crucial in velocity gauge [2(| ). The 
off-diagonal mixing term leads to to the matrix element AM = (m |(—er ■ E) [l/(E n — |r* r 7 [dE 1 /dr^)\n) + 

(m|h.c.|n), 


AM = — (p.v: 
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After the subtraction of Z-independent terms, the following two results for Ji(x) = ALL sm(2LZ)ln(\£ q + L |) 
and J 2 (x) = —d 2 Ji(x)/dx 2 supplement the analytic integrals given in Eq. 0 , 
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We can finally give the complete result for the mixing term AM, with full account of retardation, as a sum over 
virtual states |g), 
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The energy variable <? 9 is defined with respect to the reference state; i.e., if one evaluates the |m)-state admixture to 
the reference state |n), then one sets £ q = E q — E n . For excited reference states, results for both A E given in Eq. 0 
and AM in Eq. 0 contain long-range retardation tails for excited reference states, 
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where II|| and IIx are the longitudinal and transverse static polarizabilities [for the ground state, IIj_ = II|| = 11(0)]. 
The mixing term has the following long-range asymptotics, 

Asr cw 1 - 1 \ / 1 -. i\/^ ^<3 sin(2£ g Z) 3£ 2 cos(2 £ q Z) 3 £ 9 sin(2£ ? Z) 3£ q cos(2£ q Z)^ 

M = e L 0 (-^) Mmk> • tehi *l«> y- 4Z - !L -^— + 8Z3 + 16z4 — J 


,V^ A( , , . . „ 2| , (£ q cos(2£,Z) 3£ sm(2£„Z) 3 cos(2 £ q Z) 

+ e 2 ^0(-^)(m|z| g )(g|f]f-2 X 2 |n) M- 


8 Z 2 


16Z 3 


32Z 4 


( 10 ) 


•Z 5 


T ( ( m M 9) <9 Inf I 71 ) + 7 H~k) (tfkV) + | (mini k) • <flli ^l 71 ) ) + (rn|h.c.|n). 


Z> 




The results 0 and 0 will now be applied to metastable hydrogen. 

I- 

Nonretarded admixtures to metastable hydrogen .—The tential (from now on we set the elementary charge e = 1), 
results given in Eq. © and 0 have a rather involved 
analytic structure. In the short-range limit, these results 
can be compared to the static interaction of the electron 
and proton [33l , l34j with their respective mirror charges. 

This interaction leads to the following nonretarded po- 
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where we ignore terms of order 1 /Z 5 and higher [35], ]36j . 
After some tedious, but straightforward algebra, one can 
convince oneself that the terms of order Z -3 and Z -4 
are in agreement with the short-range asymptotics of the 
results given in Eqs. ([5]) and (0), i.e., in the regime Z -C 
l/£ q , which is equivalent to the limit \ —> 0. 

For close approach of the atom to the wall, the in¬ 
teraction energy is well described by the static poten¬ 
tial (ED, which necessitates a diagonalization of the 
Schrodinger Hamiltonian plus the nonretarded potential 
V (both “diagonal” interaction and Lamb shift or fine 
structure terms, as well as “mixing” terms) in the ba¬ 
sis of the 12Si/ 2 }, |2J\/ 2 ), and |2P 3 / 2 ) Schrodinger-Pauli 
wave functions with magnetic projection fj, = +1/2, to 
form the manifestly coupled states |5 i/ 2 ), |Pi/ 2 ), and 
|P 3 / 2 ). We denote the (free-space) fine-structure and the 
Lamb shift interval as 


T= 1.66 x Hr b a.u., £ = 1.61 x 1(T 7 a.u., (12) 


respectively. According to the adiabatic theorem [37H39J, 
the |<S!/ 2 ) state eigenvector has the form 


|<Si/ 2 ) ft as |25' 1 / 2 ) + ai \2P 1 / 2 ) + 0,3 |2P 3 / 2 ), (13a) 


as = 1, ai = 


15 

cz 4 


1 3 15 , nl . 

a3 = \ 1 — _ A , (13b) 


2JZ 4 


1/£»Z» l/£ 4 / 4 , 1/Z > Z > I/P 1 / 4 , (13c) 


where we ignore higher-order terms in the expansion in 
inverse powers of Z. The absolute square of the admix¬ 
ture is given by 



FIG. 1. (Color.) The modulus-squared admixtures to the 
coupled IS 1 / 2 ) state are obtained from a diagonalization of 
the potential m in the basis of |<Si/ 2 ), IP 1 / 2 ), and IP 3 / 2 ) 
states, for close approach of the atom toward the wall. The 
subscript j in Eq. ED takes on the values j = S, as well as 
j = 1/2 and j = 3/2 and denotes the state responsible for 
the admixture. As the |<Si/ 2 ) state approaches the wall, the 
initially dominant |Si/ 2 ) state contribution (solid curve, j = 
S) gradually fades and the IP 1 / 2 ) admixture (short-dashed 
curve, j = 1/2) increases, while a significant admixture of 
the IP 3 / 2 ) state (long-dashed curve, j = 3/2), is observed 
only for close approach. The atom-wall interaction energy 
becomes commensurate with the Lamb shift and fine structure 
at Z ft 84 and at Z ft 184, respectively. 


reads as 


|Pi/ 2 ) ft bs |2Si/ 2 ) + bi /2 |2Pi/ 2 ) + ^ 3/2 I 2 P 3 / 2 ), (16a) 


bs 
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6.63 x 10 15 
Z 3 
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(14) 


and bi = 1. The IP 3 / 2 ) state reads as follows, 

1 ^ 3 / 2 ) ~ c s |2Si/ 2 ) + Ci /2 |2Pi/ 2 ) + c 3 / 2 |2P 3 / 2 ), (17a) 


The one-photon decay width of the 2P state is given 
as r 2P = 6.27 x 10 s ^ = 1.51 x 10 -8 a.u., whereas 
the two-photon decay width of the 2 S state reads T 2 g = 
8.229 ^ = 1.99 x 10~ 16 a.u.. The effective decay rate 
T eff at a distance Z is 


r eff = r 2S + r 2P s = (Y 99 x 10- 16 + 1 - 01 z x 8 1 ° 8 

(15) 

We have T e ff = 2 T 2 g for Z 0 = 918 a.u.. The leading 
(nonretarded) term in the atom-wall energy shift at this 
distance amounts to —7Z,)" 3 /2 = —4.52 x 10 -9 a.u. and 
approximates both the single-particle perturbative shift 
given in Eq. (0) as well as the adiabatic energy of the 
coupled | t S 1 / 2 ) state obtained from the diagonalization of 
the potential ED to within 10%. The atom-wall energy 
at Zo is equal to —29.7 MHz and thus much smaller than 
the Lamb shift and fine structure. 

The admixture formulas for the coupled |Pi/ 2 ) state 



cs 


[3 15 

V 2PZ 4 ’ 


1 1 
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(17b) 


and of course C| = 1. For very close approach Z < 300, 
higher-order terms in the expansion of the potential 
V [see Eq. (1111) ] gradually become important. [These 
are obtained by straightforward expansion of the poten¬ 
tial ED-] Numerically determined admixtures of the cou¬ 
pled |«S' 1 / 2 ) are given in Fig. [Q they do not follow the 
asymptotic formulas for very close approach. 

Long-range tails. —The oscillatory repulsive-attractive 
dominant term in the long-range limit of the energy shift, 
for the 2 S level, goes as [see Eqs. © and ®], 


AE 2 s 


9£ 2 cos(2£Z) 
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z »2> (18) 


where we have isolated the leading term from Eq. ©, 
setting £ q = —£ and carrying the summation over the 
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virtual levels |g) = |2 Pi/ 2 ) with magnetic projections /r = 
±1/2. Somewhat surprisingly, the oscillatory terms in 
Eq. (flT)ll vanish for virtual |2 Pi/ 2 } states, so that the 
long-range coupling to the lower-lying P state vanishes. 
The leading terms in the long-range asymptotics of the 
admixture coefficients read as follows [see Eq. dai, 

3^ „ 1 . . 

0,1/2 ~ 7 tCPZ 1 2 3 4 5 ’ z ^ Z ’ ^ 19a ^ 

a 3/2 - \j\ Sin(2£ Z), Z » i . ( 19b ) 

The long-range asymptotic tail of the P 3 / 2 -state ad¬ 
mixture has an oscillatory (l/Z)-form [see Eqs. (flOl) 
and (I19bl) ]. If this tail were not suppressed by the prefac¬ 
tor C 3 7 / JF, then it could have easily provided a theoretical 
explanation for the Sokolov effect [l3l - fl6j| . because the 
(l/Z)-interaction has the required functional form to de¬ 
scribe a super-long-range term. The tail is created by vir¬ 
tual |g) = |2 Pi/ 2 ) states in Eq. (flUl) . which are energeti¬ 
cally lower than the reference 125 1 ) state. The prefactor of 
the super-long-range tail of the admixture term depends 
on details of the spectrum of the atomic system and could 
be larger for other atoms. For the Pi/ 2 -state admixture 
(term ai/ 2 ), retardation changes the 1/Z 4 asymptotics 
for short range to a 1/Z 5 asymptotics at long range. A 
full QED treatment of the admixture terms is required 
for both results recorded in Eqs. (1 19 all and (119bl) . 

Conclusions. —We can safel y co nclude that the curious 
observations reported in fl3H16l | regarding super-long- 
range 2S-2P mixing terms near metal surfaces cannot 
find an explanation in terms of a long-range effect in¬ 
volving quantum fluctuations. Both the energy shift © 
as well as the mixing term (nni) have long-range tails pro¬ 
portional to 1/Z, but the energy numerator for the 25- 
2 Pi/ 2 transition is so small (Lamb shift, a 30 cm wave¬ 


length transition) that the region in which the 1 /Z terms 
dominate is restricted to excessively large atom-wall sep¬ 
arations where the single power of Z in the denominator 
is sufficient to make the interaction energy and admix¬ 
ture terms negligible. (We should add that the inclusion 
of additional mirror charges in a cavity as opposed to a 
wall can be taken into account, in the short-range limit, 
by summing the mirror charge interactions into a gener¬ 
alized Riemann zeta function [40| and therefore cannot 
change the order-of-magnitude of the admixture terms.) 

If the observations reported in Refs. (Tiii fTT} had found 
a natural explanation in terms of a QED effect, then 
this might have had significant implications for a typ¬ 
ical atomic beam apparatus [4l| used in high-precision 
spectroscopy of atoms, potentially shifting the frequency 
of transitions involving 25 atoms in a narrow tube. For 
atom-wall separations smaller than 1000 Bohr radii, sub¬ 
stantial admixture terms are found, and the 1/Z 8 9 10 11 12 scal¬ 
ing of the effective 2 S decay rate predicted by Eq. m 
could be tested against an experiment. The clarification 
of the parity-breaking admixture terms also is impor¬ 
tant for other precision measurements in atomic physics 
which involve metastable states, such as EDM and weak- 
interaction experiments (42l - l47j . The fully retarded ex¬ 
pression for the mixing term, given in Eq. m, formu¬ 
lates higher-order QED corrections to atom-wall interac¬ 
tions beyond dipole order. Generalization of the formulas 
to, e.g., the 2 3 5i mestable state of helium is straightfor¬ 
ward. One just sums the interactions over the electron 
coordinates. 
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